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In this paper, we study the radiation from embedded sources in two-dimensional finite-size 
"photonic-quasicrystal" (PQC) slabs made of dielectric rods arranged according to a 12-fold sym- 
metric aperiodic tiling. The results from our investigation, based on rigorous full-wave simulations, 
show the possibility of achieving broadside radiation at multiple frequencies, with high-directivity 
(e.g., 15 dB) and low-sidelobes (e.g., -12 dB). We also show that leaky waves are supported by a 
PQC slab, and that the beamwidth is directly proportional to the leaky-wave attenuation constant, 
which provides a physically-incisive interpretation of the observed radiation characteristics. 

PACS numbers: 41.20.Jb, 42.25. Fx, 42.70.Qs, 71.23.Ft 



I. INTRODUCTION AND BACKGROUND 



The study of the physical properties and applications 
of aperiodically- ordered structures, related to the con- 
cept of "quasicrystals" in solid-state physics^, is ehcit- 
ing a growing attention in many fields of science and 
technology^-. In electromagnetics (EM) and optical en- 
gineering, attention has been focused on photonic qua- 
sicrystals (PQCs), mainly inspired by the theory of "ape- 
riodic tilings"Ij, i.e., arrangements of polygonal shapes 
devoid of any translational symmetry and capable of cov- 
ering a plane without overlaps or gaps. Although aperi- 
odic, these geometries can still exhibit long-range order 
and (high-order, "noncrystallographic" ) rotational sym- 
metry in a local/statistical sense^. A large body of PQC 
studies, in applications ranging from lasers to superfocus- 
ing (see Ref. jB^'for a recent review), have demonstrated 
the possibility of obtaining similar effects as those exhib- 
ited by periodic photonic crystals (PCs), with intriguing 
potentials (e.g., bandgap control, lower/multiple oper- 
ational frequencies, higher isotropy, richer defect-state 
spectrum, etc.) stemming from the additional degrees 
of freedom endowed by aperiodic order. 

One application that still remains largely unexplored 
concerns the possibility of achieving directive emission 
from embedded sources in PQCs. In periodic PCs, such 
possibility has been demonstrated in several investiga- 
tions, and its theoretical aspects are well-understood. 
In particular, directive emission has been achieved via 
the use of microcavities^ (and, possibly, slig ht lo cal 
changes in the lattice period), or planar defects^^^'^, as 
well as in defect-free configurations^ ^ by suitably de- 
signing the spatial dispersion properties. Periodic struc- 
tures interpreted as artificial materialsl^ characterized 
by an effective relative permittivity lower than one have 
also been proposed as effective devices to increase the 



directivity^JJ^^'^. Building up on the seminal work in 
Ref. 1141 the underlying phenomenology was also ex- 
plained in terms of a leaky-wave (LW) model in Refs. 
I11I12I (see also Sect. IIIB). In a similar framework, a 



detailed analysis and parameterization capable of pre- 
dicting quantitatively the directivity enhancement was 
derived in Ref. 15 for a periodic wire medium slab, in 
terms of very fast and slowly-attenuating LWs. 

For PQCs, the lack of spatial periodicity (and re- 
lated Bloch-type analytic tools) significantly complicates 
the modeling and understanding of the involved phe- 
nomenologies (see, e.g., Ref. I16ll7p . In Ref. Q8] we 
presented a preliminary study of the radiation from an 
electric line-source embedded in a two-dimensional (2-D) 
finite-size dielectric PQC slab with Penrose-type (5-fold 
symmetric, quasiperiodicP') lattice. Results showed the 
possibility of achieving moderate directivity at three fre- 
quencies. Interestingly, the lowest operational frequency 
was found to be moderately smaller than the correspond- 
ing value observed for a periodic PC slab with compara- 
ble size and number of rods, but the side- lobe level (SLL) 
was still rather high. It is worth pointing out that the 
physical interpretation of such results was not fully clear, 
although the possible applicability of LW-based models 
as in Ref. [121 was envisaged. 

In this paper, building up on the preliminary study in 
Ref. UHl we further elaborate upon this subject, by ex- 
ploring some of the open issues related to: i) the possible 
extension to other tiling geometries, ii) the physical inter- 
pretation, and in) the parametric optimization (in partic- 
ular, the SLL reduction). In this framework, we consi der 
a different tiling geometry with 12-fold symmetry^^^re- 
cently explored in connection with superfocusin^2iJ22l23j 
Via a parametric study based on full-wave numerical 
simulations^, we find several instances of directive emis- 
sion with low SLL, and address their physically-incisive 
interpretation in terms of a LW model similar to that pro- 
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FIG. 1: (a): Problem geometry, (b): Stampfli inflation rule, 
(c): Zoom of the central slab area (bullets and cross mark 
the positions of dielectric rods and line source, respectively). 
Details are explained in the text. 



posed in Ref. [TSl highlighting similarities and differences 
with the periodic PC case. 

Accordingly, the rest of the paper is organized as fol- 
lows. In Sect. ITT] we outhne the problem geometry and 
the observables utilized throughout. In Sect. |III| we 
illustrate some representative results, address their phys- 
ical interpretation, and summarize the results from our 
parametric studies. Finally, in Sect. |IV[ we provide some 
brief concluding remarks and hints for future research. 



II. GEOMETRY AND OBSERVABLES 

As in Ref. '18' we consider a 2-D (z-invariant) geometry 
(Fig. IT|a)) involving a PQC slab based on a dodecagonal 
aperiodic tiling (in the x—y plane), characterized by local 
and statistical 12-fold rotational symmetry, and gener- 
ated via the so-called Stampfli inflation rulesi^ illustrated 
in Fig. flTb). The construction, to be iterated up to the 
desired tiling extension, starts from a dodecagon made 
of square and equilater-triangle tiles (thick solid line), 
which is first scaled-up by an infiation factor \/3+2 (thick 
dashed line), and subsequently substituted at each ver- 
tex of the inflated version (see Refs. I20I21I22I23I24I25I26I 
for similar and different 12-fold symmetric PQCs). The 
PQC slab is then generated by cutting a rectangular por- 
tion of size L X h from a suitably large tiling, and placing 
circular dielectric rods (in vacuum) at the tile vertices; no 
defects are assumed. The rods have relative permittivity 
er = 12 and radius r = 0.138a, with a being the lattice 
constant (tile sidelength). A time-harmonic (exp(— iwt)) 
electric line-source is placed in vacuum nearby the slab 
center (see the zoom in Fig. Illc)). 

As meaningful observables, we consider the normalized 



"local density of states" (LDOS^ 

p(rs,cj) =4:lm{G{rs,rs;uj)}, 



as well as the directivity 
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the normalized intensity 
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and the SLL 
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(1) 
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(4) 



with G denoting the PQC Green's function relating the 
current to the magnetic vector potential, E^-^ and E^q 
the far-flelds radiated by the line-source in the presence 
and absence, respectively, of the PQC, and B the main 
lobe angular region. The above observables, equal to 
one in the absence of the PQC, are computed via a 
well-established full-wave numerical technique based on 
a multipolar Fourier-Bessel expansionN. In order to facil- 
itate comparisons with previously published results, we 
also consider the 3dB beamwidth A03ds, deflned as the 
angular region where the radiation is higher than -3dB 
of the maximum, i.e. (for a symmetric radiation pattern 
with main lobe at cj) = 90°), 

\E{f{7T/2 + A,P3dB/2,Lu)\ = \E{f{7r/2,u;)\/V2. (5) 

In the 2-D scalar scenario of interest here, the LDOS 
admits an intuitive interpretation as the normalized to- 
tal power (per unit-length along the z direction) radi- 
ated by the line-source^, and represents a meaningful 
bandgap indicator, whereas the directivity and the nor- 
malized intensity are indicative of the angular conflne- 
ment and power density enhancement (as compared to 
vacuum) of the radiated far-fleld. Attention is focused 
on achieving broadside (0 = 90°) directivity, although 
the structure also radiates into the y < half-space. 



III. DIRECTIVE EMISSION FROM A 
DODECAGONAL PQC SLAB 

A. Representative Results 

Figures [2ja)-[2][c) show the observables in ([l}-([3]) as 
a function of the normalized frequency a/Xo (with Ao 
denoting the vacuum wavelength). From the LDOS re- 
sponse (Fig. p[a)), one can identify a main bandgap (rep- 
resented by low values of p) and several minor dips. Po- 
tentially interesting conflgurations are characterized by 
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FIG. 2: Normalized LDOS (a), broadside directivity (b), and 
broadside normalized intensity (c), as a function of the nor- 
malized frequency a/Xo, for a dodecagonal PQC slab (L = 
97.8a, h = 4a, 435 rods, cf. Fig. flj;), with source at i = 0, 
y — —0.5a. 




FIG. 3; Radiation patterns (in the upper half-space) for a 
configuration with the same parameters as in Fig. [2l at three 
normalized frequencies (marked as vertical dashed lines in Fig. 
[2|. (a), (b), (c): a/Ao = 0.239,0.317, and 0.736, respectively. 



high- amplitude peak coincidences in the broadside direc- 
tivity (Fig. l2Fb)) and normalized intensity (Fig. |2jc)). 
These coincidences correspond to situations where the 
peaks in the radiated power (Fig. l2Fc)) are hkely accom- 
panied by radiation patterns with narrow beams, and 
thus high directivity, in the (p = 90° direction. The verti- 
cal dashed lines mark three such representative configura- 
tions, whose radiation patterns (normalized |i?/-'^(0, a;)p 
vs. 0) are displayed in Figs. l3Fa)-l3jc). In all three 
cases, one should observe the narrow peak at broadside 
accompanied by sidelobes of lower level. It should be em- 
phasized that, as clearly visible from the sharp peaks in 
Fig. [21 the above configurations are intrinsically narrow- 
band. At the lowest frequency (a/Ao — 0.239), the direc- 
tivity is about 9dB (with 3dB beamwidth Ai/is^jb = 4°), 
and the SLL (Fig. Isla)) is moderately high (A w —6 
dB), comparable with the values observed in Ref. [18] 
for a Penrose- type PQC. However, at the higher fre- 
quencies a/Ao = 0.317, 0.736, the radiation patterns 
look considerably cleaner, with directivities of 13.7dB 
{A(j)3dB = 2.35°) and 10.4dB (A03dB = 3.44°), respec- 
tively, and the SLL reaches values ^ — lOdB that may 
be acceptable in many practical applications. As com- 
pared with the Penrose tiling geometry in Ref. [TSJ the 
dodecagonal PQC here turns out to provide comparable 
operating frequencies and lower SLL. 

It is worth noticing that direct comparisons of above 
results with the periodic PC cases in Refs. 1415161718151 
should be addressed cautiously, in view of the differences 
in the underlying radiation mechanisms and related para- 
metric ranges. In this connection, we highlight that the 
configurations in Refs. I4I5I6I7I8W1 are based on either 
defected or defect-free thick structures, which essentially 
exploit the spatial filtering capabilities of PCs around 
band-edge frequencies in order to achieve directive emis- 
sion. Conversely, our configuration is based on a defect- 
free thin-slab geometry, and, as one can observe from Fig. 
[2] the operational frequencies are not necessarily close to 
the band edges. As a reference for qualitative compar- 
isons, we observe from Ref. [5|that the radiation pattern 
produced by a line-source embedded in a defect-free pe- 
riodic PC thick-slab made of 24 x 16 alumina (e^ = 9.61) 
rods of radius r = 0.136a (in the experiment, it was cho- 
sen a = 1.1cm and r = 1.5mm), at a frequency of 13.21 
GHz (i.e., a/Ao = 0.484), would exhibit a 3dB beamwidth 
A03ds = 8°, and SLL values ~ — 20dB. Therefore, from 
a qualitative viewpoint (and bearing in mind the afore- 
mentioned parametric differences), the performance ex- 
hibited by our proposed configuration (at the higher fre- 
quencies) turns out to be comparable, with moderately 
smaller beamwidth (higher directivity) and higher SLL. 

We also would like to point out, that in general, the 
directivity of thick structures is governed by geometrical- 
optics ultra-refraction phenomena, as explained in Ref. 
1131 for a metamaterial with a rectangular lattice and in 
Ref. [28] for its equivalent homogenized low-permittivity 
material. On the other hand, thin slabs are able to pro- 
vide higher directivitiy because the physical mechanism 
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FIG. 4: (Color online) Internal field intensity maps pertaining 
to the radiation patterns in Fig. Is] (due to symmetry, only the 
a; > slab-half is shown), (a), (b), (c): a/Ao = 0.239,0.317, 
and 0.736, respectively. 
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FIG. 5: Field-magnitude cut at y = 2.03a for the configura- 
tion in Fig. W[c) (corresponding to the radiation pattern in 
Fig. jSTc)). Note the semi-log scale and the linear fit (shown 
dashed), which highlight the exponential field decay along x. 



is slightly different, as explained in Ref. i28, There, it 
is shown that, because of the small thickness, a (sin- 
gle) dominant LW can be excited. In Ref. '^S' and in 
Sect. III-D of Ref. 29l the two mechanisms (geometrical- 
optics-based for thick structures, and LW-based for thin 
structures) have been quantitatively compared, showing 
that the latter is capable of yielding a higher directivity. 



I11I12I131 which can be effectively parameterized via a LW 
model (see Refs. ll4)15|29|30t3T] for more details). For 
slabs made of periodic structures, it has been observed^ 
that the field excited by a source and evaluated at a cer- 
tain distance from the vacuum-slab interface is mainly 
represented by a LW mode contribution. Such distance 
from the interface is required in order to observe a homo- 
geneous field distribution that is not strictly related to 
the local structure details described by evanescent waves 
(associated with high-order Bloch harmonics) that decay 
away from the interface. Along the interface, at a fixed y, 
a general description of the field in terms of algebraically- 
decaying spatial waves and guided (bounded and leaky) 
modes is given in Ref. [HI Sect. Ill of Ref. [HI and 
Sect. V of Ref. [301 In certain structures (see, e.g., 
those shown in Ref. I5CT)| . the spatial wave is the domi- 
nant contribution. However, in other cases (which can be 
properly designecP^, the spatial wave is comparatively 
weakly excited, and the near-field Es at an equivalent 
aperture close to the vacuum-PC interface is therefore 
dominated by a LW mode with complex propagation con- 
stant k^w = P + ia, i.e.. 



Es{x, y) ^ E„{y) exp {if3\x\) exp (-a|a;|) , (6) 



where /3 is the phase propagation constant, and the atten- 
uation constant a describes the exponential decay along 
X. If the attenuation constant a is much smaller than 
the free space wavenumber fep = 2n/Xo, the field at the 
vacuum-PC interface is sustained for a long distance, ex- 
tending for several wavelengths. If also the phase prop- 
agation constant /? is much smaller than the free space 
wavenumber (i.e., /3 <C kg), we derive that the field at the 
interface exhibits an essentially equi-phase profile, which 
is an important aspect for high-directivity radiators. It 
is important to note that close to the vacuum-PC inter- 
face (at distances smaller or comparable with the lattice 
constant) the field oscillates around the average value 
represented by ^ because of the effects produced by the 
nearby scatterers (for periodic structures, these oscilla- 
tions are represented by high-order Bloch harmonics). 



B. Physical Interpretation 

In this Section, we show that the radiation mechanism 
underlying the observed results is related to the excita- 
tion of a LW that travels along the PQC slab. A LW is a 
guided wave that radiates continuously across a guiding 
structure while propagating along it; the power leakage 
results in an exponential attenuation of the wave as it 
propagates along the structure, even in the absence of 
losses. 

The radiation mechanism of our PQC thin slab turns 
out to be somehow similar to those illustrated in Refs. 



As anticipated, another important outcome of this 
study is the evidence that LWs can be supported by a 
PQC slab as well. In order to gain such evidence, we 
show that the field at the vacuum-PQC interface can be 
described by the simplified model in (|6]) . As a first exam- 
ple. Fig. |4]shows the intensity field-maps inside the PQC 
slab at the three representative frequencies, from which 
the presence of a modal structure is clearly visible. 

The LW propagation constant k^w is very useful to pa- 
rameterize the characteristics of the radiated field, which 
can be readily computed via Fourier transform of the LW 
aperture field in ^. We can therefore write the LW ra- 
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FIG. 6: (Color online) Comparison between the actual radi- 
ation pattern at normalized frequency a/Xo — 0.736 (blue 
solid curve, cf. Fig. [3[c)) and LW-model prediction for 
(3 = a = 0.099/a (red dashed curve, cf. (tsl). 
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where A is a normalization constant. From (|7|, it can 
been shown that a radiation peak at broadside is ob- 
tained wherP^'291 p <^ ko and (3 < a. Moreover, the 
highest directivity at broadside is fou nd when the opti- 
mum condition /3 « a is satisfiecP^ESl^ jj^ which case the 
radiation pattern becomes 
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From this parameterized expression of the radiation pat- 
tern, we can derive the 3dB beamwidth. 



A03ds ~ 2 arcsin 



\/2a 
fco 



2^/2c 
fco 



(9) 



At variance of periodic PCs and artificial materials, in 
our aperiodic geometry, the LW propagation constant 
khw cannot be systematically computed from a disper- 
sion equation. Nevertheless, we can show that the ob- 
served radiation features are consistent with the predic- 
tions of a simplified LW model as in (l6|-(|9]). To this aim, 
for the operational frequency a/Ao — 0.736, Fig. [sjshows 
a cut of the field magnitude at an aperture y — 2.03a 
close to the vacuum-PQC interface. Thanks to the semi- 
log scale and the dashed linear fit in the plot, an expo- 
nential decay is clearly observable, with irregular oscilla- 
tions around the linear fit ascribable to the inherent local 
structure of the PQC slab interface. The LW attenua- 
tion constant a can readily be estimated from the slope 
of the log-scale linear fit, yielding a ~ 0.099/a, which is 
consistent (a ^ fco) with the observed high directivitjffSl 
Due to the aforementioned lack of a dispersion equation, 
we are not able to compute (3. Nevertheless, the ob- 
servation from Fig. ^c) that at the chosen frequency 



FIG. 7: As in Fig. [3|a), but with source at a; = 0, y = 1.55a. 



a/Ao = 0.736 the directivity is actually maximum sug- 
gests that the condition (3 ~ a is satisfied (see the dis- 
cussion after ( [7|) . This is confirmed by the comparison, 
shown in Fig. |6J between the actual radiation pattern 
(cf. Fig. l3{c)) and the LW prediction from (|8|, which 
evidences a good agreement in the angular region near 
the maximum; the two curves are almost coincident over 
a 5 dB dynamic range. From a more quantitative view- 
point, Eq. (pi estimates the corresponding beamwidth 
^4>3dB = 3.46°, in excellent agreement with the actual 
value of 3.44°. We can thus conclude that the observed 
radiation properties are indeed consistent with the pre- 
dictions of the simplified LW model in (|6|-((9]), and hence 
a LW with the appropriate characteristics (namely, with 
amplitude and phase profiles varying on a much larger 
length scale than the vacuum wavelength) is mainly re- 
sponsible for shaping the main beam of the pattern. 

Another interesting feature offered by the LW param- 
eterization is the information about the transverse size 
(along the x direction) that is necessary to produce a 
required beam radiation. Once the beamwidth A(j)3dB 
is specified, the required a — (3 attenuation constant is 
then found via ([9|. From ([6|, one can then determine 
how large the structure needs to be for having a negligi- 
ble field at its edges, so that truncation effects would not 
produce pattern distortion in the main lobe region. 



C. Summary of Parametric Studies 

From the above physical interpretation, one gathers 
that directive low-SLL radiation can be achieved by ex- 
citing a proper LW mode in the PQC slab. In order to 
gain further insight into the role of the various param- 
eters involved in the excitation process, we carried out 
a comprehensive set of parametric studies, whose main 
results can be summarized as follows. 

• The symmetric-cut configuration in Fig. [l] turns 
out to be the best suited for exciting a LW radi- 
ating at broadside. Different cuts from the orig- 
inal dodecagonal tiling, where either the x— or 
y— symmetry is broken, can be utilized for pointing 
the main beam in different directions (as already 
observed in Ref. .181). 



• Increasing the slab thickness, we also observed vari- 
ations in the radiation characteristics, though with- 
out the clearly-pred ictable be havior observable in 
the periodic PC caseEMIIMIl. in the periodic PC 
case, it is known that increasing too much the slab 
thickness, other LW modes may appear and dete- 
riorate the radiation pattern. In this connection, it 
should also be noticed that, unlike in periodic PCs, 
thickness changes in PQC slabs intrinsically involve 
termination changes, and so the two effects (thick- 
ness and termination) are not easily discernible. 

• The rods radius value (r = 0.138a) utilized 
throughout was arrived at via parameter scanning. 
Slight departures from this value (e.g., r = 0.15a) 
were observed to yield moderate deterioration in 
the SLL, up to the complete disappearance of the 
directive response at the higher frequencies. 

• The source position was found to be a rather ef- 
fective parameter for adjusting the LW excitation, 
so as to optimize the radiation pattern. As an ex- 
ample, the radiation pattern in Fig. p[a) was op- 
timized by moving the source vertically (so as to 
preserve the a;-symmetry and, hence, the broad- 
side radiation) via parameter scanning. Figure u\ 
shows the optimized radiation pattern, which rep- 
resents the best tradeoff response {D = 15.5dB, 
S = 14.8dB, and A = -12dB) observed for a do- 
decagonal PQC in our study. 



IV. CONCLUSIONS 



In this paper, we have presented a study of the radi- 
ation properties of line-sources embedded in 2-D finite- 
size defect-free dielectric dodecagonal PQC slabs. Via a 
comprehensive full-wave numerical study, we have iden- 
tified certain parametric configurations which give rise 
to directive low-SLL broadside radiation. These re- 
sults extend to a different PQC geometry our previous 
observations^'^, confirming the possibility of obtaining di- 
rective low-SLL radiation at multiple frequencies. More- 
over, we have shown that also for the case of aperiodically- 
ordered PQCs the observed radiation characteristics are 
consistent with the predictions of a LW-based simplified 
model, and have explored the role of the key geometrical 
parameters. 

The above results, and their LW-based physical inter- 
pretation, pave the way for a systematic parameteriza- 
tion and design of directive radiators based on PQC slabs. 
In this framework, current and future investigations are 
aimed at a deeper understanding of the modal structure 
supported by PQC slabs, and the refinement of the LW- 
based modeling, so as to come up with more quantitative 
design rules. Also of interest is the experimental verifi- 
cation of the phenomenon, as well as further studies of 
different configurations (e.g., air-hole-type) and PQC ge- 
ometries. 
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